Influence of plural scattering on the quantitative determination of spin and orbital 
moments in electron magnetic circular dichroism measurements 
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Recent quantitative measurements of the orbital to spin magnetic moment ratio ruL/ms in elec- 
tron magnetic circular dichroism (EMCD) experiments have given a rriL/rns ratio that is larger than 
commonly accepted values. We demonstrate here that plural scattering may noticeably influence 
the niL/ms ratio. An equation is derived which describes its influence as a function of the spectral 
integrals of the plasmon scattering region and zero-loss peak. The influence of the electron-plasmon 
scattering can be removed when electron energy-loss spectra of the ionization edge are deconvoluted 
by the low-loss signal. For a bcc-Fe sample we obtain mL/ms= 0.04 after plasmon removal. We 
conclude that the plural scattering should be considered when extracting quantitative information 
from EMCD measurements. 
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I Electron magnetic circular dichroism (EMCD) is a 
young experimental techniquCfi which in principle al- 
lows to extract atom-specific spin and orbital magnetic 
moments from electron energy-loss spectra by applying 
EMCD sum rulesj^i^ This technique, as compared to its 
counterpart x-ray magnetic circular dichroism (XMCD)^, 
brings a promise of nanometer or even sub-nanometer 

I spatial resolution^ with bulk sensitivity using a stan- 
dard transmission electron microscope (TEM). 

Recently a few works reported the first EMCD evalu- 
ations of the orbital to spin moment ratios 771^/7715 — 

I for bcc-Fe 0.09 ± 0.03^ 0.08 ± O.Oljl and for hcp-Co 
0.14 ± 0.03^). It should be noted that all these values 
overestimate the 777^/7775 ratios as measured by Chen et 

• al^ using XMCD (0.043 for bcc-Fe and 0.095 for hcp-Co). 

' In this paper we demonstrate that a possible explanation 
of this overestimation is based on plural scattering, i.e., 
low-loss excitations accompanying the core-level ioniza- 
tion processes. 

In our recent worki we reported measurements of 
energy-resolved diffraction patterns on a 20 nm bcc-Fe 
, thin film, which were acquired in two-beam (2BC) and 
three-beam (BBC) case orientations of the sample. Based 
on symmetry argument o^'^° the SBC is preferred and a 
new double difference procedure has been introduced for 
extraction of the magnetic component of the signal. The 
electron energy-loss spectra (EELS) were in Ref. pro- 
cessed for every pixel in reciprocal space and thus evalu- 
ated by extracting data from maps of the signal. Instead 
of working with 777 l /itis values from individual pixels it 
is possible to use averaged spectra from the same set of 
dataii, which should lead to the same results under sym- 
metrical conditionsiS. Here, spectra were obtained from 
the data cube of Ref. (after application of point blem- 
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FIG. 1: Illustration of the influence of plural scattering. Top: 
convolution of the low-loss region, consisting of zero-loss peak 
and plasmon peak (at -Epl), with a peak at El^ leads to a 
double-peak structure. Bottom: the same procedure applied 
on Z/2 peak. If E^^ is close to El^ + -BpL, then the L2 peak 
is enhanced by the plural scattering accompanying the L3 
excitation. 

ish removal, cross-correlation and rotation, see Refs. [10 
for details) by selecting two virtual apertures of 70x70 
pixels (corresponding to a square of size G(200)/2), ac- 
cording to the optimum position and size determined in 
Ref. 0. Each spectrum was thereafter background sub- 
tracted, normalized in the post-edge region and fitted us- 
ing the same models as described in Ref. @. The method 
of averaged spectra improves the signal to noise ratio and 
is therefore preferable here for studying the influence of 
plural scattering. 

It is assumed that the effect of the plural scattering 
can be modelled by a convolution of a net (single scat- 
tering) near-edge spectrum with the spectrum from the 
low-loss regioufi^ as illustrated in Fig. [TJ An L3 peak at 
energy-loss E^^ is convoluted with the low-loss spectrum, 
which has a double-peak structure: an intense zero-loss 
peak (ZLP) and a relatively weak plasmon peak at Epj^. 



Result of this convolution is again a double-peak struc- 
ture with main peak at Ej^^ and an additional peak at 
Eph + El^ . The latter appears due to double excitation 
events, i.e., when one fast electron excites a plasmon and 
a L3 core-level transition. The energy-integral of this ad- 
ditional peak is proportional to the intensity of the L3 
peak and to the ratio k of energy-integrals of plasmon 
peak (Api) and ZLP (Azlp), k = Apl/Azlp- It is im- 
portant to observe that if this peak overlaps with the L2 
edge, then the L2 edge is enhanced by low-loss inelastic 
scattering events accompanying the L3 excitation. 

The low-loss spectra, measured from the bcc-Fe sample 
at four different thicknesses, arc shown in Fig. [S] Note 
that the peak-to-pcak ratio between the plasmon peak 
and the ZLP is only around 1%. However, the spectral 
integral of the plasmon region (within the energy-loss 
range 10 to 35 eV) is approximately 10% and 12% of the 
ZLP area at thicknesses t/X = 0.18 and 0.23, respectively 
[t is thickness and A is the mean free pathii). Note 
that the EpL is approximately at 22 eV with full- width 
at half-maximum above 20 eV, Fig. [2l and the energy 
difference of Lj, and L2 edges in Fe is cca 13 eV. Therefore 
in bcc-Fc the plasmon excitations accompanying the L3 
excitation indeed overlap with the L2 peak. 

A simplified estimation of the influence of plural scat- 
tering on the Tahjms ratio can be obtained in the fol- 
lowing way: the convolution of near-edge spectra with 
low-loss spectra enhances the L2 peak proportionally to 
the ratio fc and the intensity of the L3 edge.— Assuming 
that symmetry conditions for detector positions required 
by the EMCD sum-rules^ are fulfilled, the spectra at de- 
tector positions H- and - are schematically given by: 
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L3 L2 net plural scattering 

W+mI + \N^^^^2+ 'k{Ni+Mij] (1) 
Ni-Mi + [N2 - Ah + k{Ni - Ml)] (2) 

^3 L2 total 



where N, M denote the non-magnetic and magnetic com- 
ponent of the signal, indices 1 and 2 correspond to lower 
and higher energy-loss edges, respectively. For late tran- 
sition metals Mi and M2 have typically opposite signs. 
The EMCD difference spectrum is, consequently: 
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2Mi + [2M2 + 2kM{\. 



(3) 



The effect of plural scattering does not cancel out when 
taking the difference spectrum because of the different 
intensity of the magnetic terms for the two detector po- 
sitions. 

When we insert the energy integrated values for indi- 
vidual edges into the sum- rule expression for the m l l^ng 
ratiof^ we obtain after some algebra the relation between 
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FIG. 2: Low-loss EELS spectra of bcc-Fe measured at sample 
thicknesses t/A = 0.18, 0.23, 0.32 and 0.52. At i/A = 0.32 the 
beam passes also through a V/Fe superlattice and at tjX = 
0.52 additionally through MgO substrate. The zero-loss peak 
height was normalized to one and centered at eV for all 
thicknesses. 



'clean' TO^/ms value, m, and the 'convoluted' one, u: 

{2L-l)u + L{L-l)k/3 + kLu 

u — 

2L-l-kL-3kuL/{L-l) 

(2L - l)u - L(L - l)k/3 - kLu 

2L-l + kL + 3k{iL/{L - 1) 



(4) 
(5) 



where L orbital quantum number of final states {L ~ 2 
for p ^ d and L = 3 for d f). At u = the plural 
scattering has no influence on mi/ms ratio, i.e., u ~ u. 
It can be shown that the value u = corresponds 
to Ml = 0. This explains the /c-indcpendcncc of u, see 
Eq. Q. 

By expanding into first order in k and omitting a term 
proportional to the square of the clean rriL/ms ratio 
(negligible for late transition metals) we obtain simpler 
equations for £2,3 edges: 
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The second term in Eq. ^ predicts that if A: = 10%, 
then the m^/ms ratio obtained by straightforward ap- 
plication of sum rules is enhanced by 0.022. The fac- 
tor in parentheses has a smaller effect under the as- 
sumptions made here, giving an enhancement by ca. 
13%. Applying Eq. ^ to the value obtained in Ref. 
(?7iL/ms)conv. = 0.08 leads to {mL/ms)c\ca.n = 0.05, 
which is a considerable change. 

Figure [3] shows the prediction of Eq. © for differ- 
ent strengths of the plasmon region, fc, thus revealing 
the relation between convoluted and deconvolutcd (clean) 
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FIG. 3: (color online) a) Qualitative effect of plural scattering 
for different values of the area ratio fc between the plasmon 
region and ZLP, Eqs. (j4| and ([5]). b) mi/ms value as a func- 
tion of a; = M2/M1 (ratio of the integrated dichroic signals 
from L2 and L3 edges). 



rriL/ms ratio. From the raw spectra, convoluted values of 
the areas are available, namely Mi and M2 = M2 + kMi. 
When also k is known, from Fig. [S] 'clean' rriL/ms ratios 
can be estimated. According to the EMCD sum rules^^ 
the clean or convoluted rriL/ms ratio is a function of 
X = M2/M1 or M2/M1, respectively. For the ^2,3 edges 
this gives m^/ms = |(1 + x)/{l — 2x), which is shown 
in Fig. [1). 

To study the effect on a quantitative level, wc per- 
formed a plasmon removal procedure on our experi- 
mental data taken in the SBC orientation. Fourier- 
ratio deconvolutioni^ii^iii^ was performed in Digital 
Micrograph"'"'^, using a zero-loss modifier as reconvolu- 
tion method. To test the sensitivity of the extraction 
of the data we positioned the set of boxes first at the 
Thales circle positionsji and subsequently at the optimal 
positions with maximum relative dichroic signal^^ Each 
set comprises of four spectra, one in each quadrant of 
the diffraction plane. It should be noted that the ZLP 
could not be acquired under exactly the same illumina- 
tion conditions as the near-edge spectra due to very high 
intensity of the electron beam in the telefocus mode^, 
which could damage the CCD camera. 

In Fig. 0^ we show the averaged spectra from all four 
quadrants at the optimal positions. Optimum position- 
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FIG. 4: (color online) Influence of the plasmon removal on 
energy-loss and dichroic spectrum, see text for details, a) 
Averaged normalized spectra from all four quadrants (see leg- 
end); b) the same spectra after plasmon removal; c) The dif- 
ference of spectra from panels a) and b). In all panels: the 
red curve with hatching relates to the double difference signal 
(DD) and the grey curve to its energy-integral. Values p and 
q enter the sum rules expression of Ref. 0. 



ing of the box in combination with the double difference 
procedure results in a strong EMCD signal, reaching al- 
most 30% of the L3 edge peak height. Analyzingthis sig- 
nal by means of peak- fitting as described in Ref. (a we ob- 
tain (mi/TO5)conv. = 0.09, in agreement with the value of 
Ref. [Tl where the value was extracted from ttil /ms ratio 
maps. Figure |4)d shows the same spectra after plasmon 
removal using the low energy-loss spectrum measured at 
t/X = 0.23, which corresponds to the sample thickness at 
which the data were obtained. The resulting change is 
relatively weak, as could be expected taking into account 
the small plasmon area as compared to the ZLP. How- 
ever, the higher energy-loss tails end at an intensity of 
approximately 0.85 instead of unity as in the normalized 
spectra of Fig. [4^. 

Figure \^ displays the effect of plasmon removal as a 
difference between spectra from panel a) and panel b). 
The onset of the enhancement starts approximately at 
10 eV above the onset of the edge, in agreement with 
the low energy-loss spectra, Fig. [5J where the 0-10 eV 
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Thales circle Optimum pos. 

Orientation raw data deconv. raw data deconv. 

Up-down, right 007 Om 007 002 

Up-down, left 0.10 0.06 0.10 0.05 

Double difference 0.09 0.04 0.09 0.04 



TABLE I: Summary of the mi/ms values obtained in SBC 
geometry from raw spectra and from spectra after plasmon 
removal, both at the optimal placement and Thales circle po- 
sitions, see text for details. 

region is dominated by the ZLP. At about 723 eV there 
is a shoulder which originates from the same feature of 
the plasmon peak at around 15 eV. The broad humps 
around 735 eV and 750 eV are derived from L3 and L2 
edges smeared by the wide plasmon peak. These features 
are, compared to £2,3 peak positions, shifted up in energy 
by approximately 20-25 eV. This shift is associated to the 
position of the plasmon peak, around 22 eV above zero 
energy-loss. 

Since the plasmon removal will alter also the double- 
step background^ which in turn can influence the L3 
edge area, we applied the peak fitting procedure de- 
scribed in Refs. 100 to deconvolved spectra. Fig. IHd. The 
results are summarized in the Table HI There we show 
nT-L/fTT-s values in BBC orientation, using a square inte- 
gration box of O.5G200 X 0.56*200 Bit optimal placement 
and Thales circle positions, prior to and after the plas- 
mon removal. The up-down difference for left and right 
side refers to the use of a horizontal mirror axisi^ The 



double difference combines both horizontal and vertical 
mirror axesi^ It effectively cancels out deviations in the 
TfiL/ms ratio, induced by asymmetry and misalignments 
in the experiment j^ii^ After plasmon removal we obtained 
for bcc-Fe a mi/m^ ratio of 0.04 ± 0.01 in this setup. 

We checked the method on the only other EMCD 
rriL/ms value available in the literature, the hcp-Co^. 
We estimate that for a 18 nm thick hcp-Co sample, the 
value for k (in the same energy range 10-35 eV) would be 

0.1. Applying the Eq.[5]on value {mL/ms)conv = 0.14 
gives {rni^/mg) clean ~ 0.10, which agrees with XMCD 
value 0.095^ within its error bar. 

We have demonstrated that the plural scattering has a 
non-negligible influence on quantitative information ex- 
tracted from EMCD spectra. Therefore, a measurement 
of the low energy-loss spectra and subsequent plasmon 
removal from the near-edge spectra is a necessary part of 
the quantitative EMCD analysis. A qualitative formula 
for estimation of the error introduced by plural scatter- 
ing was derived. The correction is found to scale ap- 
proximately linearly with the ratio of spectral integrals 
of the plasmon region and ZLP. By applying the plas- 
mon removal procedure to spectra measured in Ref. [T] we 
obtain tol/tos = 0.04 ± 0.01, which is in good agree- 
ment with values 0.043 and 0.044 obtained by XMCD^ 
and gyromagnetic ratioi^, or 0.062 by neutron scattering 
experiments^^ respectively. 
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The mean free path A for 300 keV electrons in Fe is roughly 
around 100 nm.— The pure ZLP was obtained by us- 
ing different peak fitting models available in the Digital 
Micrograph^'^ software, such as refiected tail or combina- 
tion of Gaussian/Lorentzian (resulting in a FWHM of ap- 
proximately 2 eV). 

More precisely one should take the energy-integral of the 
plasmon peak within the suitable energy range, over which 
the near-edge spectra are peak-fitted, here energy-loss of 
ca. 10-35 eV. 



